Capillary adhesion between elastic solids with randomly rough surfaces 
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I study how the contact area and the work of adhesion, between two elastic solids with randomly 
rough surfaces, depend on the relative humidity. The surfaces are assumed to be hydrophilic, and 
capillary bridges form at the interface between the solids. For elastically hard solids with relative 
smooth surfaces, the area of real contact and therefore also the sliding friction, are maximal when 
there is just enough liquid to fill out the interfacial space between the solids, which typically occurs 
for (Ik ~ 3/irms, where cJk is the height of the capillary bridge and /irms the root-mean-square 
roughness of the (combined) surface roughness profile. For elastically soft solids, the area of real 
contact is maximal for very low humidity (i.e., small dx), where the capillary bridges are able to 
pull the solids into nearly complete contact. In both case, the work of adhesion is maximal (and 
equal to 27COS&, where 7 is the liquid surface tension and 9 the hquid-solid contact angle) when 
(Ik » hrms, corresponding to high relative humidity. 



1. Introduction 

When two solids are in close contact capillary bridges 
may form at the interface, either as a result of liquid-like 
contamination layers (e.g., organic contamination from 
the normal atmosphere) or water condensation in a hu- 
mid atmosphere, or intensionally added thin fluid layers, 
e.g., thin lubrication films. For wetting liquids strong 
negative pressure will prevail in (short) capillary bridges, 
which will act as an effective relative long-ranged attrac- 
tion between the solids [ij. In many cases, in particu- 
lar for hard solids and with rough surfaces, the contri- 
bution to the wall-wall attraction from capillary bridges 
may be much larger than the contribution from the di- 
rect wall- wall interaction, e.g., the van der Waals inter- 
action between the solids. The capillary bridge mediated 
wall-wall interaction has a huge number of important ap- 
plications, e.g^ for granular materials^, insect or tree 
frog adhesion[3, 0, head/disk systems @ and micro- 
electromechanical systems (MEMS) , where they may 
trigger permanent adhesion and device failure. 

The influence of capillary bridges on adhesion is well 
known to all of us. Thus, it is possible to build sand 
castles from humid or slightly wet sand but not from dry 
sand or sand flooded with water. Similarly, it is well 
known that very flat surfaces, such as those of gauge 
blocks (steel blocks with the roughness amplitude of or- 
der ~ 25 nm when measured over a macroscopic area, 
e.g., 50 cm^), adhere with strong forces resulting from 
capillary bridges formed from organic contamination or 
water. Thus it may easily be shown that the force to 
separate two gauge blocks may strongly increase if one 
exposes the surfaces to (humid) breath (for strong ad- 
hesion, the fluid layer (organic or water) should be at 
least several times larger than the surface root-mean- 
square roughness amplitude, i.e in the present case of 
order 100 nm, or more; see below). Finally, the fact that 
a fly can walk on a vertical glass window is due to cap- 
illary bridges formed at the tip of many thin hair-like 
fibers, which cover the attachment organs of the fly, see 
Fig. [TJ Without the fluid (injected via channels in the 
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FIG. 1: The attachment organ of a fiy. At the tip of the 
hairs on the attachment pads occurs thin elastic plates. The 
fly inject a liquid to the space between the plate and the 
substrate. The resulting capillary bridge makes it possible 
for the fly to adhere to surfaces with large roughness. To 
very rough surfaces, the fly can also "adhere" using the claws 
on the legs. Adapted from S. Gorb, with permission. 



fibers) adhesion would probably be impossible to most 
surfaces because too much elastic energy is necessary to 
bend the tip of the fibers into atomic contact with the 
rough substrate, which is necessary for strong adhesion 
without fluids. [The tips of the hair covering the attach- 
ment pads of some lizards (e.g., Geckos) and most spiders 
are much thinner that those of a fly, and can easily be 
bent to make atomic contact even to very rough surfaces; 
in these cases no liquid is injected by the animal into the 
contact area.] 

Capillary adhesion between solids with randomly 
rough surfaces has been studied so far mainly using the 
Greenwood- Williamson contact mechanics theory's']. In 



this theory the asperities on the rough surfaces are ap- 
proximated by spherical cups and the long-range elastic 
coupling is neglected. It has recently been shownjl, 
that for surfaces with roughness on many length scales, 
the GW theory (and other asperity contact theories such 
as the theory of Bush et al[ll|) fail qualitatively to de- 
scribe e.g., the area of real contact and the interfacial 
separation as a function of the load. 

When two elastic solids with rough surfaces are 
squeezed together, the solids will in general not make con- 
tact everywhere in the apparent contact area, but only 
at a distribution of asperity contact spots 12, 3, 
The separation m(x) between the surfaces will vary in 
a nearly random way with the lateral coordinates x = 
(x, y) in the apparent contact area. When the applied 
squeezing pressure increases, the average surface separa- 
tion u = (u(x)) will decrease, but in most situations it 
is not possible to squeeze the solids into perfect contact 
corresponding to m = 0. In thermal equilibrium, capil- 
lary bridges are formed at the wall-wall interface in those 
regions where the separation w(x) is below the Kelvin 
length (iK, which depends on the relative humidity and 
on the liquid contact angles with the solid walls. 

Here I will present a general theory of capillary ad- 
hesion between elastic solids with randomly rough sur- 
faces. The study is based on a recently developed the- 
ory for the (average) surface separation -0 as a func- 
tion of the squeezing pressure p. The theory shows that 
for randomly rough surfaces at low squeezing pressures 
p ~ exp(— au//irms) where a ~ 2 depends (weakly) on 
the nature of the surface roughness but is independent 
of p, in good agreement with experiments [3]. The GW 
contact mechanics theory (and the more accurate theory 
of Bush et al) instead predicts p ~ u~°'exTp[—b{u/ h^^s)'^] 
(where a and b are positive numbers) in str ong dis- 
agreement with numerical simulations [itI [isl . Il9j and 
experiment [3l ■ 

The theory presented below is based on solid and fluid 
continuum mechanics. Thus it cannot be strictly applied 
when the Kelvin distance becomes of order the molec- 
ular size, i.e., smaller than 1 nm. The limiting case of 
molecular thin fluid films must be treated using atom- 
istic methods, e.g.. Molecular Dynamics (see, e.g., Ref. 

MEH). 



2. Theory 

Consider the frictionless contact between two elastic 
solids with the Young's elastic modulus Ei and E2 and 
the Poisson ratios i^i and 1^2- Assume that the solid sur- 
faces have the height profiles /ii(x) and /i2(x), respec- 
tively. The elastic contact mechanics for the solids is 
equivalent to those of a rigid substrate with the height 
profile /i(x) = /ii(x) -I- /i2(x) and a second elastic solid 
with a flat surface and with the Young's modulus E and 
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FIG. 2: Capillary bridges formed at two asperity contact ar- 
eas. At thermal equilibrium, the radius of curvature of the 
capillary bridge tk is given by the Kelvin radius. The capil- 
lary bridges will exert an attractive force _Fa on the block. The 
sum of the capillary force _Fa and the external load Fq = poAo 
(where Aq is the nominal surface area) must equal the repul- 
sive force arising from the area of real contact between the 
solids. 




FIG. 3: An rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification it appears 
as if complete contact occur, but when the magnification is 
increased it is observed that in reality only partial contact 



the Poisson ratio ly chosen so that [22 
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Consider an elastic block with a flat surface in contact 
with a hard rough substrate. We consider humid condi- 
tion (vapor pressure Pv) and assume that the liquid wet 
the surfaces. In this case in the asperity contact regions, 
liquid capillary bridges will form (see Fig. [2]), where the 
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FIG. 4: An asperity contact region observed at the magnifica- 
tion ^. It appears that complete contact occur in the asperity 
contact region, but upon increasing the magnification to the 
highest resolution (magnification ^i) it is observed that the 
solids are separated by the (average) distance u(C)- 



meniscus radius rx is given by the Kelvin equation 



fcBrln(Pv/Psat)' 



(2) 



where vq is the molecular volume in the liquid. Thus, the 
thickness of the liquid film is given by 



dK — ?'k(cos6'i + cos6'2), 



(3) 



where 9i and 02 are the liquid contact angles on the two 
solid surfaces. 

For water at T = 300 K, 7 = 0.073 J/m^ and vq ~ 
3 X 10-29 m^ so that 



0.53 nm 

"ln(Pv/Psat 



Thus, if water wet the surfaces (i.e., 9i = 62 = 0) the 
thickness of the liquid film is given by the Kelvin length 



1.06 nm 

'ln(Pv/Psat)' 



(4) 



We now use the contact mechanics formalism devel- 
oped elsewhere 17, 2^ 24, 25, [2^, where the system is 
studied at different magnifications see Fig. |3l When 
the system is studied at the magnification ^ it appears 
as if the contact area (projected on the xy-plane) equals 
A{(), but when the magnification increases it is observed 
that the contact is incomplete, and the surfaces in the 
apparent contact area A{(^) are in fact separated by the 
average distance u{(), see Fig. HI Let mi(C) be the (aver- 
age) height separating the surfaces which appear to come 
into contact when the magnification decreases from C. to 
C — AC, where is a small (infinitesimal) change in the 



magnification. ui{() is a monotonically decreasing func- 
tion of and can be calculated from u{Q and A{C) using 
(seeRef. ^) 



(5) 



We assume that liquid occurs in the apparent contact 
areas when the separation wi(C) is smaller than (or equal 
to) the Kelvin distance dx- Assume that this occurs for 
the magnification ( — so that 



ui(Ck,Po) 



(6) 



where we have indicated that the separation ui depends 
on the nominal pressure pq. In the liquid bridge is a 
negative pressure p — ~pk with pK = '^j/d-K- The liquid 
occupy the (projected) surface area A^ — A{(k) — A{(i), 
where ^(Ci) is the area of real contact observed at the 
highest magnification (i . Thus the attractive force 

F,=PKAA^^[AiCK)-AiCi)]. (7) 

We define p^ = Fa/Ao. We will calculate A{Q by using 
a mean-field type of approximation, where instead of in- 
cluding the non-uniform distribution of capillary forces 
acting at the block-substrate interface, we assume that 
the effective squeezing pressure is p = Po + Pa, where 
Po = Fq/Aq is the applied squeezing pressure (which is 
negative during pull-off) (see Fig. [5] and also Sec. 5). 

Let us apply a pull-off force Fq = — Fpuii to the block. 
We define the work of adhesion per unit area as 



w = — / du Ppull(w) = / du [pa 
^0 JuQ Juo 



p{u)]. (8) 



p{u) is the repulsive pressure from the substrate at the 
separation u — u between the average substrate surface 
plane and the average position of the bottom surface of 
the block, uq is the equilibrium separation when -Fpuii = 
0, i.e., p{u(i) = Pa- 
Let us now study the limiting case when the space 
between the solids is filled with hquid, i.e., no dry area. 
We also assume that the area of real contact ^(Ci) is 
negligible compared to the nominal contact area ^o- In 
this case the attractive pressure Pa = 27/c?k must be 
balanced by the repulsive asperity contact pressure which 
for separation u > ft-rms is given by 



p{u) — f3E*e 



(9) 



where a and (3 are numbers which depend on the surface 
roughness but which are independent of p and of the 
elastic properties of the solids 18]. Thus 
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a 27 



(10) 



(11) 



Since the capillary pressure is Pa = 27/dK for u < 
and zero otherwise, we get work of adhesion 



du [pa - p{u)] 



= 3— (ctK - Mo) 
«K 
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Using (11) this gives 



w — 2j 



1 



-log 



27 



For self-affine fractal surfaces we have[18j a 
(3 « 0.4go/irms giving 



w « 27 



' ^rms T 
1 - -T^log 



2dK 



27 



(12) 



(13) 



2 and 



(14) 



where go is the roll-off (or cut-off) wavevector of the sur- 
face roughness power spectrum. 

The analysis above has assumed that thermal equilib- 
rium occurs at any moment in time during pull-off, which 
requires a static or slowly propagating debonding crack. 
In this case liquid will condense or evaporate at the in- 
terface in such a way as to always maintain the Kelvin 
radius tk = (^k/Z for the radius of curvature of the liquid 
meniscus. However, during fast pull-off negligible con- 
densation or evaporation occurs, and the fluid volume at 
the interface, rather than the Kelvin radius, will be con- 
stant during pull-off. If A{u) is the area covered by liquid 
when the average surface separation is u, then volume 
conservation requires AqUq = A{u)u or A{u)/Aq = uq/u. 
In this case 



du Pa{u) = 

Uq J Uo 



, 27 A(u) „ 
du — = 27. 

u Aq 



When thermal equilibrium occurs during pull-off the 
same integral becomes 27(1 — uo/^k)- Thus, when 
mq/'^k << 1, the whole interface is filled with fluid, and 
the adiabatic pull-off (constant meniscus radius) and fast 
pull-off (constant fluid volume) give nearly the same re- 
sult for the work of adhesion, assuming that one may 
neglect viscous energy dissipation in the fluid during the 
fast pull-off (which is responsible for suction-cup type of 
effective adhesion). 

In general, for most solids with covalent, ionic or metal- 
lic bonds E « 10^"'^ Pa, and in the normal atmosphere 
(where dx is of nanometer size) capillary adhesion will 
only be observed if /irms is at most a few nanometer. In 
fact, assuming that (14) is valid we get w = if 



, J ^rms, / qohimsE*dc 

d^<d.^ —log ^ 



(15) 



Since dc only depends logarithmically on the parameters 
go, E* and 7, the critical Kelvin distance (or humidity). 



rms = 6 




logioq (1/m) 
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FIG. 5: The surface roughness power spectrum of two sur- 
faces. The surface with the root-mean-square (rms) roughness 
6 fim is self affine fractal with the fractal dimension Dt = 2.2. 
The other surface has the rms roughness 2.4 nm and is as- 
sumed to be randomly rough. 



below which the adhesion is very small, depends weakly 
on these parameters. For example, if ft-rms ~ 2 nm, and 
with the roU-of wavevector go ^ 10'' m~^ (as for the lower 
curve in Fig. [5]) we get log(go/irms£'*c'K/27) « 5 so that 
dc ~ 5 nm corresponding to 80% relative humidity. 

As another application of (14), let us study the case 
where E* = 1 GPa (as typical for glassy polymers) in 
contact with a hard rough substrate with /ir,ns = 1 Mui 
and go = iC m~^. Using 7 « 0.1 J/m^ (14) gives that 
w > only if cZk > 10 /im. More generally, (14) shows 
that capillary adhesion will only manifest itself as long 
as the the Kelvin length c?k is (at least) a few times the 
root-mean-square amplitude of the (combined) substrate 
roughness profile. This statement only holds for clasti- 
cally hard enough solids-for very soft solids the negative 
capillary pressure can pull the solids into close contact in 
which case (14) is no longer valid, see Sec. 3.2. 

3. Numerical results 

We will apply the theory to two cases, namely the con- 
tact between elastically hard solids, e.g., silicon as rele- 
vant for microelectromechanical systems (MEMSs), and 
elastically soft systems such as rubber or biological adhe- 
sive pads. For the case of hard solids we will assume that 
the (combined) rough surface has the surface roughness 
power spectrum C(g) shown in Fig. [5l bottom curve. 
This power spectrum was obtained from the surface to- 
pography, ft.(x), measuredj27j for a polysilicon surface, 
using [28j: 



Ciq) 



1 



(2n) 



d^x {h[^)h{0))e-''^-'- 



The root-mean-square (rms) roughness of this surface is 
2.4 nm. For the case of elastically soft solids we will as- 
sume that the combined surface has the power spectrum 
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FIG. 6: The stress as a function of the logarithm of the (av- 
erage) separation u during separation. The area under the 
curves determines the work of adhesion. For a hard surface 
with the root mean square roughness /irms = 2.4 nm in contact 
with an elastic solid (with the Young's modulus E — 82 GPa 
and Poisson ratio v — 0.22) with a flat surface. Results are 
shown for the relative humidities (RH) 0.95, 0.925, 0.9, 0.875. 
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FIG. 7: The work of adhesion as a function of the inverse of 
the height dx of the capillary bridges. For a hard surface with 
the root mean square roughness /irms ~ 2.4 nm in contact with 
an elastic solid (with the Young's modulus E — 82 GPa and 
Poisson ratio u — 0.22) with a flat surface. The square data 
points are the calculated results for the relative humidities 
(RH) 0.95, 0.925, 0.9, 0.875, and the solid line a fit to the data. 
The line denoted by "approx" is given by (13) with a — 1.87 
as calculated [itI from the measured surface roughness power 
spectrum. 

shown in Fig. [5l top curve. This is the power spectrum of 
a self affine fractal surface with the rms roughness 6 /im 
and the fractal dimension D[ = 2.2. 

3.1 Elastically hard solids 

In Fig. [HI we show the stress as a function of the (aver- 
age) distance between the surfaces, u, during separation 
for several different relative humidities (RH), 0.95, 0.925, 
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FIG. 8: The work of adhesion as a function of the relative 
humidity. The circle are experimental data (from Ref. f29D 
and the squares calculated results. In the calculation we have 
assumed a hard surface with the root mean square rough- 
ness /irnis = 2.4 nm in contact with an elastic solid (with the 
Young's modulus E = 82 GPa and Poisson ratio u = 0.22) 
with a flat surface. The square data points are the calculated 
resuhs for the relative humidities (RH) 0.95, 0.925, 0.9, 0.875. 
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FIG. 9: The work of adhesion as a function of the relative 
humidity. In the calculation we have assumed hard surfaces 
with the root mean square roughness /irms = 2.4 and 3 nm 
in contact with an elastic solid (with the Young's modulus 
E — 82 GPa and Poisson ratio f = 0.22) with a flat surface. 
The square data points are the calculated results and the lines 
smoothing cubic splines. 

0.9 and 0.875. The elastic solid has the Young's modulus 
E = 82 GPa and the Poisson ratio v — 0.22. The area 
under the curves in Fig. [S] determines the work of adhe- 
sion which is shown in Fig. [T) The numerical results in 
Fig. [7] are in relative good agreement with the (approxi- 
mate) analytical result (14) (curve denoted by "approx" 
in Fig. El). 

In Fig. [5] we compare the work of adhesion, as a func- 
tion of the relative humidity, with the experimental data 
of DelRio et al. [1^ [s^l , obtained from microcanterlever 
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FIG. 10: The area of real contact as a function of the Kelvin 
capillary width (Ik- For a hard surface with the root mean 
square roughness /irms = 2.4 nm in contact with elastic solids 
with the Poisson ratio v — 0.22 and the Young modulus E = 
82 GPa. The (nominal) squeezing pressure p = 4 MPa. 
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FIG. 11: The dry and wet area as a function of the Kelvin 
capillary width (Ik- For a hard surface with the root mean 
square roughness hims = 2.4 nm in contact with elastic solids 
with the Poisson ratio v — 0.22 and the Young modulus E = 
82 GPa. The (nominal) squeezing pressure p = 4 MPa. 



experiments. The circles are experimental data and the 
squares calculated results (from Fig. [7]). The experimen- 
tal system is for polysilicon surfaces with the combined 
surface having the root-mean-square (rms) roughness (29j 

2.3 nm, while in the calculation we have assumed the 
same system as above where the rough surface has the 
rms roughness /irms ~ 2.4 nm. 

In Fig. [HI I show the calculated work of adhesion as 
a function of the relative humidity. In the calculation 1 
have assumed hard surfaces with the root mean square 
roughness hrms = 2.4 (from Fig. [8|) and 3 nm. The 
power spectrum of the second surface was obtained from 
the original power spectrum (bottom curve in Fig. [5]) 
by scaling it with a factor of (3/2.4)^. The slope of the 



rms = 3 nm line is about 25% larger than for the rms = 
2.4 nm curve. This is just the ratio between the two rms 
values and agree with the prediction of the (simplified) 
theory, Eq. (14), which shows that the slope of the work 
of adhesion curve scales approximately linearly with the 
rms-value. This prediction is very different from the GW 
theory prediction which predict much larger change in 
the work of adhesion. 

In the present case the theory predicts that the work 
of adhesion per unit area, w, vanishes for RH < 0.87 (see 
Fig. [5|) or, equivalently, for the capillary heights < 
8 nm (see Fig. [?])• In the experiment w is indeed very 
small for RH < 0.87 but not zero (see Fig. [8]), and we 
attribute this difference between theory and experiment 
to two different effects: 

1) Finite size effect: The theory is for an infinite system 
while the experimental system is finite, and in the present 
case in fact quite small (the crack tip process zone has 



a total area of only ~ 100 /J.m^, see Ref. |3ll)- A finite 
pull-off force due to capillary bridges will in reality always 
occur, even in the limit of very low RH where (for a small 
system and solids which high elastic modulus) a single as- 
perity (or just a few asperities) may be in contact-in this 
case a (small) capillary bridge can form at the asperity 
giving a non-zero pull-off force. For a finite-sized system 
this will result in a non-zero work per unit area, w, to 
separate the solids (which may be non-negligible for mi- 
cro and nano-sized systems), but in the thermodynamic 
limit (infinite sized system, as assumed in the theory), 
w would of course vanish. The small contact between 
small, elastically hard, solids also imply that there will 
be large fiuctuations in the pull-off force between differ- 
ent realizations of the same system. This has indeed been 
observed by Zwol et al[32|, who found that the pull-off 
force (at low relative humidity) varied by a factor of ~ 2 
(or more) , when the same microsized object was brought 
into contact with the same substrate surface at different 
locations (see also Sec. 5). 

2) Van der Waals interaction: One can easily show 
that the long-ranged van der Waals interaction will al- 
ways gives a non-zero work of adhesion. The reason for 
this is that the (repulsive) contribution to the wall-wall 
interaction from the elastic deformation of asperities de- 
cay as ~ exp(— a-u/Ziriiis) (or faster, if finite size effects 
are taken into account) with increasing wall separation 
M, while the (attractive) van der Waals interaction decay 
slower as ~ 1 /u^ and the total interaction will always be 
attractive for large enough separation u. With only the 
capillary bridges the wall-wall interaction is predicted to 
be repulsive (for all wall- wall separations) at low rela- 
tive humidity. Combining the van der Waals interaction 
and the capillary contribution will give a non-zero work 
of adhesion even for small relative humidity, which still 
depend on the relative humidity. I plan to study this in 
more detail in the future. 

It is interesting to compare (14) with the predictions 



of the Greenwood-Williamson (GW) theory of contact 
mechanics. As pointed out before, this theory fails qual- 
itatively when roughness occurs on many length scales, 
and DelRio et al. found that the GW theory strongly 
underestimate the work of adhesion. The main reason 
for this failure is not the asperity approximation (which, 
however, also is a severe approximation in the present 
case) but the neglect of elastic coupling between the con- 
tact regions. 

The separation of two solids usually occurs via interfa- 
cial crack propagation, which depends on the work of 
adhesion. However, the sliding friction is determined 
mainly by the area of real contact. In Fig. [TO]we show the 
area of real contact as a function of the Kelvin length c?k- 
The (nominal) squeezing pressure is p = 4 MPa. Note 
that the area of real contact is maximal at ~ 7 nm, 
which correspond to the point where there is just enough 
fluid to fill the space between the solids. This is illus- 
trated in Fig. [Tl] which shows the wet (and dry) area 
as a function of the Kelvin length cZk- Note that for 
c^K > 7 nm the interface if filled with liquid. In this 
case the negative pressure 27/dK will prevail (nearly) ev- 
erywhere at the interface and the area of real contact 
will therefore be proportional to l/dx, which is indeed 
the dependence of A/Aq on cZk observed in Fig. [TUl for 
cIk > 7 nm. In fact, we can accurately describe the be- 
havior of the contact area for large and small cIk- Thus, 
when no liquid occur at the interface, and since 

the applied (squeezing) pressure po is very small the area 
of real contact is linearly related to the po according to 

A 

-r = apo 

where a can be calculated from the surface roughness 
power spectra and the effective elastic modulus as de- 
scribed elsewhere [illii]. We get a = 3.1 x 10"^° Pa"^ 
which gives the limit A/Aq ~ 1.24 x 10^"^ as cIk 0, in 
good agreement with Fig. [TUl 

In the opposite limit of large (nearly) the whole 
interface is filled by fluid so that the effective pressure 
acting on the block is the sum of the applied pressure po 
and the capillary pressure 27/dK giving 

This relation (denoted "fully wet interface" ) is shown in 
Fig. [TU] and agree very well with the numerical result for 
c^K > 7 nm. 

3.2 Elastically soft solids 

Consider the contact between an elastically soft solid 
with a flat surface and a hard, randomly rough, substrate 
with the power spectrum shown in Fig. \E\ upper curve. 
This system exhibit more complex and interesting behav- 
ior than the case of elastically hard solids studied in the 
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FIG. 12: The area of real contact as a function of the loga- 
rithm of the Kelvin capillary width For a hard surface with 
the root mean square roughness ftrms = 6 /im in contact with 
elastic solids (with a flat surfaces) with the Poisson ratio 0.5 
and several different Young modulus, E = 3, 10, 30, 100 and 
300 MPa. The (nominal) squeezing pressure p = 0.1 MPa. 
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FIG. 13: The area of real contact as a function of the loga- 
rithm of the average width, dav, of the water layer between 
the surfaces. For a hard surface with the root mean square 
roughness /irms ~ 6 fim in contact with elastic solids with the 
Poisson ratio 0.5 and several different Young modulus, _E = 3, 
10, 30, 100 and 300 MPa. The (nominal) squeezing pressure 
p = 0.1 MPa. 

last section. Thus, we find that when the relative humid- 
ity decreases, resulting in a higher negative pressure in 
the capillary bridges, there may be a strong increase in 
the area of real contact because of elastic deformations 
of the solid walls, which are pulled into closer contact by 
the capillary adhesive forces. This is in sharp contrast 
to the case of hard solids (with relative smooth surfaces) 
studied above, where the contact area was maximal when 
there was just enough liquid to fill the space between the 
(almost) undeformed solid walls. 
Fig. [T2] shows the area of real contact as a tunc- 




FIG. 14: The dry and wet area as a function of the loga- 
rithm of the average width, dav, of the water layer between 
the surfaces. For a hard surface with the root mean square 
roughness /irms ~ 6 pm in contact with elastic solids with the 
Poisson ratio 0.5 and the Young modulus E — 10 MPa. The 
(nominal) squeezing pressure p — 0.1 MPa. 
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FIG. 15: The stress as a function of the (average) separation 
u during separation. The area under the curves (for positive 
pull-off stress) determines the work of adhesion. For a hard 
surface with the root mean square roughness /irms = 6 fj,m 
in contact with an elastic solids (with the Young's modulus 
S = 3, 10 and 100 MPa and Poisson ratio = 0.5) with a flat 
surfaces. Results are shown for the Kelvin distance = 1 M^i 
(corresponding to the relative humidity RH ~ 0.999). 



tion of the logarithm of the Kelvin capillary width cIk, 
for a hard surface with the root mean square roughness 
^rms = 6 /im, in contact with elastic solids (with flat 
surfaces) with the Poisson ratio 0.5, and several differ- 
ent Young modulus, E = 3, 10, 30, 100 and 300 MPa. 
The (nominal) squeezing pressure p = 0.1 MPa. Note 
that for low enough RH (i.e., for small enough capillary 
distance (Ik) the area of real contact increases rapidly 
with decreasing RH. This result from the strong capillary 
bridges pulling the solids into close contact. As expected, 
this effect start at higher and higher RH as the elastic 



modulus of the solids decreases. For the stiffest solid, 
E = 300 MPa, the area of real contact does not increase 
when the RH is reduced even to the point that the height 
of the fluid bridges is of molecular size dx ~ 0.5 nm (be- 
low which the continuum theory fail). 

Fig. [13] shows again the area of real contact but now 
as a function of the logarithm of the average width, dav, 
of the water layer between the surfaces. This figure illus- 
trates an interesting effect: as the RH decreases, in some 
range of RH the amount of liquid between the surfaces 
increases as the RH decreases, i.e., water from the at- 
mosphere condenses at the interfacial region between the 
solids. The reason for this is that when the RH decreases, 
the surfaces are pulled together which imply that the wet 
interfacial area [where 'u(x) < c?k] increases, and this in- 
crease is so fast that the volume of liquid between the 
surfaces increases in spite of the fact that the capillary 
height decreases. 

This effect is also illustrated in Fig. [T3] which show 
the wet (and dry) area as a function of the logarithm of 
the average width, dav, of the water layer between the 
surfaces. 

In Fig. [15] we show the stress as a function of the 
(average) interfacial separation u, during pull-off. The 
area under the curves (for positive pull-off stress) deter- 
mines the work of adhesion. Results are shown for elastic 
solids with a flat surfaces, and with the Young's modulus 
E = S, 10 and 100 MPa and the Poisson ratio v ^ 0.5. 
The results are for the Kelvin distance d}^ — 1 /im [cor- 
responding to the relative humidity (RH) « 0.999]. Note 
that for the most stiff solid E = 100 MPa the interaction 
is purely repulsive, i.e., the work of adhesion vanish. For 
the elastically softer solids the work of adhesion is finite, 
but the pull-off stresses are much smaller than for the 
hard smooth surfaces studied in Sec. 3.1., while the dis- 
tance over which the pull-off force is non-negligible are 
much longer, both effects reflecting the much larger fluid 
film thickness (or Kelvin distance d^) in the present case. 
The work of adhesion (i.e., the area under the curves in 
Fig. [13 for the E = 3 MPa case is w = 0.11 J/m^ and 
for the E ^10 MPa case we get w = 0.062 J/m^. 

4. Applications 

We present a few applications of the results presented 
above. Consider first rubber friction on a wet glass sur- 
face in the context of wiper blades. After rain or car 
wash, the driver of a car does not immediately stop the 
motion of the wiper blades. In this case, as the water is 
removed by wiping and evaporation, as a function of time 
a high friction peak, with a friction considerably higher 
than the dry one, may be observed. This time period of 
enhanced friction is denoted as the tacky regime, and may 
prevail for several seconds. Experiment [33| have shown 
that in the tacky regime the attraction from water capil- 
lary bridges between the rubber and the glass substrate 



pull the rubber into close contact with the substrate so 
that the area of real contact is even larger in the tacky 
regime than for the perfectly dry contact. Since the ap- 
plied pressure in the rubber-glass nominal contact area 
is very high in wiper blade applications, typically of or- 
der MPa, the additional contribution from the capillary 
bridges must be very large, of order MPa, in order to 
explain the strong increase (typically by a factor of 2) 
in the friction. This implies that the height of the capil- 
lary bridges d < 100 nm. This also implies that the wiper 
blade-glass rubber friction may be considerably enhanced 
at humid condition, where capillary bridges are formed 
spontaneously. 

The hight friction in the tacky region can block a wiper 
system. Capillary adhesion can be reduced by surface 
treatment of the rubber. Thus, halogenation result in 
a rubber surface layer which is elastically much stiffer 
than the bulk (by up to a factor of 100, see Ref. [s^ D 
and will strongly reduce the friction coefficient in the 
tacky region, e.g., from 3 (no halogenation) to 1 (at 15% 
halogeneous concentration in the surface region) [35| . The 
surfaces of wiper blades are almost always exposed to 
halogenation. However, after long time of use this layer 
of modified rubber may be removed by wear. Thus, the 
friction coefficient for very worn-out blades may exceed 
3. 

Capillary bridges may also give rise to strongly en- 
hanced friction for lubricated rubber applications, if the 
lubrication film is very thin everywhere (e.g., even at the 
boundary of the nominal contact area). Thus, in a re- 
cent experiment [3EI, a rubber block was slid on a smooth 
steel surface lubricated by a drop of oil. In this case the 
oil film where everywhere very thin d w 300 nm, and 
the rubber-steel surface friction was observed to increase 
from ~ 1 to ~ 3 when the oil drop was added to the metal 
surface. In this experiment the applied nominal pressure 
was rather small, about 0.01 MPa. The capillary adhe- 
sive force can be estimated using 2j/d « 0.3 MPa where 
we used 7 « 0.05 J/m^. This is much larger than the 
applied nominal squeezing pressure, but because of the 
finite sliding speed (0.1 m/s) one cannot assume that all 
the fluid get completely squeezed out from the asperity 
contact regions [37|, and therefore the increase in the slid- 
ing friction will be smaller than indicated by the increase 
in the effective normal load. 

Capillary adhesion is very important in biological ad- 
hesive systems used for locomotion. This has been stud- 
ied in details for tree frogs [1, [3| and stick insects [11]. 
These animals use smooth adhesive pads which are built 
from non-compact materials and are elastically very soft. 
To adhere to surfaces the animals inject a wetting liq- 
uid into the contact area. Experiments for stick insects 
have shown that the nominal frictional shear stress for 
sliding against smooth glass surfaces increases when the 
fluid fllm thickness decreases, e.g., during sliding long 
distances. During repeated sliding on the same surface 



area, the friction decreased continuously as more liquid 
accumulated on the surface - typically the frictional shear 
stress dropped from 0.15 MPa to 0.05 MPa as the film 
thickness increased. Nevertheless, as expected, the nom- 
inal pad-substrate surface contact area did not change 
much since the work of adhesion (which determines the 
contribution to the contact area from the adhesional in- 
teraction) may be nearly constant (equal to 27) if enough 
liquid occurs at the interface, as indeed expected in the 
present case (note: the insects usually need to move on 
surfaces much rougher than the glass surface and must 
therefore inject much more liquid at the pad-substrate 
contact area than is necessary for the smooth glass sur- 
face). The pull-off force from the glass surface was nev- 
ertheless observed to increase when the film thickness 
decreased during repeated contact and pull-off. This is 
most likely a viscosity effect: during separation of two 
closely spaced solids separation by a thin fiuid layer, 
fiuid must fiow towards the center of the contact area, 
and during "fast" pull-off this generates a strong nega- 
tive pressure (effective adhesion) in the fluid because of 
viscous dissipation in the fluid. In this case, the thinner 
the fluid fllm, the larger peak force will be, as seen in the 
experiments. 

As a flnal application of the theory presented above, let 
us consider the contact between the toe pad of a tree frog 
and a smooth glass plate. This system has been studied 
in great detail in Ref. 4] . At the toe pad-substrate inter- 
face occurs a wetting liquid, which is likely to be of order 
~ 1 fim thick at the edges of the toe pad-substrate con- 
tact regions (see Ref. [5|). Thus, the negative capillary 
pressure 2j/d « 0.1 MPa. Because of the low effective 
elastic modulus of the toe pad, it deforms elastically so 
that a large fraction (about 50%) of the toe pad comes 
into close contact (less than 5 nm) with the glass surface, 
in spite of the fact that the natural height fluctuations 
of the toe pad surface (on the length scale of ~ 10 ^m) 
is of order 1 /xm. The area of real contact manifests 
itself experimentally as a finite static friction force. 

5. Discussion 

The theory presented in Sec. 2 is a mean-field type 
of theory, where the adhesive force from the capil- 
lary bridges is taken into account by adding an adhesive 
pressure Pa = -Fa/^o to the external squeezing pressure 
Pq = Fq/^o acting on the block {Fq is the normal load, 
which is assumed to act uniformly on the top surface of 
the block). This is likely to be an accurate approximation 
as long as a large fraction of the non-contact interfacial 
region is filled with liquid, e.g., for high relative humid- 
ity. However, for very low relative humidity and for large 
surface roughness, liquid bridges may only occur close to 
the outer edges of the asperity contact regions. If the 
regions occupied by fluid are very small compared to the 
diameter of the contact and non-contact regions, then we 



can consider the fluid filled regions as crack-tip process 
zones, and this limit can be studied using the adhesion 
theory developed elsewhere [i^. [3§|. In this theory the 
interface is studied at different magnifications, and an ef- 
fective interfacial binding energy (per unit area), 7efi(C), 
is introduced. If one assumes that the continuum the- 
ory of capillary forces is still valid for the small capillary 
bridges which occur at low RH then 7eff(Ci) — 2^cosd. 
However, most likely the continuum theory is not valid 
for narrow (molecular sized) capillary bridges, and one 
may have to treat 7eff(Ci) = ^7 ^'^ empirical param- 
eter to be determined directly from experimental data. 

In the theory developed in Sec. 2 we always assumed 
thermal equilibrium. However, it is known that the for- 
mation of capillary bridges in a humid atmosphere is a 
thermally activated process, with a continuum of activa- 
tion energies and hence relaxation times. Thus, in gen- 
eral it may take very long time to reach (or come close 
tojthermal equilibrium. In the study by Maarten et 
al|3l| for microcantilevers, at the relative humidity 0.3, 
the average crack length (defined as the non-contact part 
of the cantilever) decreased with increasing times from 

700 yum initially to « 550 fim after 2.5 h, to ~ 400 
after 5 h, and to « 375 after 10 h. After this time no 
further decrease of the crack length was observed. In gen- 
eral, the increase in adhesion and friction with increasing 
time, as a result of thermally activated formation of cap- 
illary bridges, has been studied in detail, and is of great 
importance in many systems of fundamental or applied 
interest 0, EH, El]. 

Another limitation of the theory presented in Sec. 2 is 
that it is only valid if most of the relevant repulsive wall- 
wall interaction occurs for separations u that are small 
enough so that, for the actual physical system, there 
are enough asperity contact regions (to obtain enough 
self-averaging) that the analytical theory can be applied. 
That is, the analytical contact mechanics theory is for an 
infinitely system (thermodynamic limit), and for a ran- 
domly rough surface (with a Gaussian probability height 
distribution) there will always be infinite high asperi- 
ties, and contact between two solids will occur for any 
(average) separation u between the solids. In fact, for 
large u we have the exact result p ~ exp(— au/Zimis), 
which shows that a repulsive pressure p will act between 
the solids for arbitrary large separation u. However, for 
finite-sized systems, the highest asperities have a finite 
height, which for macroscopic systems t ypi cally is of or- 
der ~ 10/irms (see Appendix A in Ref. [23|), but in the 
context of MEMS may be much smaller due to their small 
physical size. In fact, de Boer has estimated the high- 
est asperities in the crack tip process zone to be only 
w 3.7/irms- It is clear that the analytical theory is only 
valid if the important repulsive contribution to the work 
of adhesion occurs for separations u < 3.7/irms- In our 
applications to microcantilever adhesion most of the rel- 
evant repulsive interaction occurs for u < 3.2/ii-ms so the 



theory is likely to be at least semiquantitatively correct. 

In the GW theory a similar problem related to the 
height of the highest asperities, but in addition the exact 
form of the tail of the height distribution for large separa- 
tion matters a lot for the contact mechanics and in partic- 
ular for the relation between u and p. However, this is an 
artifact of the GW theory, and when the long-range elas- 
tic coupling is taken into account as in the theory of Pers- 
son, the contact mechanics becomes much less dependent 
on the details of the height distribution. This has been 
verified by Finite Element Method calculations [l^ for 
polymer surfaces [3] where, in spite of the fact that the 
height probability is rather non-Gaussian for large as- 
perity heights (see Fig. 14 in Ref. [i3|)j the contact 
mechanics obeys the usual behavior with the area of real 
contact proportional to the load and the average inter- 
facial separation depending on the squeezing pressure as 
p ~ exp(— au/ferms), in good numerical agreement with 
the Persson theory [l7|. This can be explained by the 
fact that because of long range elastic deformation, not 
only the asperities close to the top of the Gaussian distri- 
bution will make contact, but a larger range of asperity 
heights will be involved (which is the reason for why the 
the asymptotic relation p ^ exp(— aw//irms) is exponen- 
tial rather than Gaussian (reflecting a Gaussian height 
distribution) p ~ exp[— 6(u//irms)^] as in the GW the- 
ory). This is illustrated in Fig. [16] for a situation where 
this fact is particularly clear. 

The analytical theory presented in Sec. 2 for the limit- 
ing case of an interface completely filled with fiuid (e.g., 
high relative humidity) indicates that the most impor- 
tant property of the rough surface is the root-mean- 
square roughness /iims, while other properties such as 
(for self affine fractal surfaces) the fractal dimension D{ 
and the upper and lower cut off wavevectors qi and qo 
are much less important because the parameter a in the 
relation p exp(— au/ft-rms) is very insensitive to these 
quantities This is likely to be the case also for non- 
fractal surfaces. This fact is consistent with experiments 
of van Zwol et al[32|, who observed that surfaces with 
very different qo but the same root-mean-square rough- 
ness /iims ~ 1-5 nm, exhibited the same capillary ad- 
hesion. We conclude that the rms roughness gives the 
predominant effect on the adhesive force due to capillary 
bridges. 

In the theory presented in Sec. 2 we have neglected 
plastic deformation of the solids and the disjoining pres- 
sure due to adsorbed (water) layers. The problem of the 
influence of plastic yielding on the contact mechanics was 
studied in Ref. ^29J where it was found to have a very 
small influence on the capillary adhesion. I have analyzed 
the same problem using the more accurate contact me- 
chanics theory I have developed. I find that the surface 
roughness with wavelength longer than ~ 0.2 /xm under- 
goes negligible plastic deformation, while some plastic de- 
formation occurs for roughness at shorter length scales. 
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FIG. 16: (a) In the Greenwood-Williamson (GW) contact 
mechanics model only local (asperity) elastic deformations are 
included, (b) In reality the elastic deformation is long ranged, 
and asperities which never could be in contact in the GW 
theory because they are too low, could be in contact when the 
elastic deformations is fully included in the analysis. For this 
reason the GW theory is much more sensitive to the exact 
form of the probability distribution of asperity heights for 
large asperity height, than in a more accurate analysis which 
include the long-range elastic deformations. 

Plastic deformation will tend to smoothen the surfaces 
and hence reduce the stored elastic energy. This in turn 
will enhance the capillary adhesion. I will report on this 
study elsewhere. 

The presence of adsorbed water layers on the hy- 
drophilic polysilicon surfaces used in the experiment [29| 
will also play a role in the capillary adhesion process. 
Unfortunately, a water adsorption isotherm does not cur- 
rently exist for polysilicon surfaces, which makes it im- 
possible to definitively determine the role of adsorbed wa- 
ter layers on the capillary adhesion problem. However, 
the fact that we obtained good agreement with experi- 
ment neglecting this effect for high relative humidity, in- 
dicates that it may not be very important (for RH > 0.9) 
for the microcantilever applications discussed above. 

6. Summary and conclusion 

I present a general theory for how the contact area and 
the work of adhesion between two elastic solids with ran- 
domly rough surfaces depends on the relative humidity. 
The surfaces are assumed to be hydrophilic, and capil- 
lary bridges form at the interface between the solids. For 
elastically hard solids with relative smooth surfaces, the 
area of real contact and therefore also the sliding friction, 
are maximal when there is just enough liquid to fill out 
the interfacial space between the solids, which typically 
occurs for /ik ~ S/irms, where /ik is the height of the cap- 



illary bridge and /irms the root-mean-square roughness of 
the (combined) surface roughness profile. For elastically 
soft solids, the area of real contact is maximal for very 
low humidity where the capillary bridges are able to pull 
the solids into nearly complete contact. In both case, 
the work of adhesion is maximal (and equal to 27COS6', 
where 7 is the liquid surface tension and 6 the liquid-solid 
contact angle) when >> /irms, corresponding to high 
relative humidity. 

The theory is compared to experimental data for mi- 
crocantilever structures. The theory is in good agreement 
with the data while the classical Greenwood- Williamson 
theory fail qualitatively. I also present applications to 
rubber wiper blades, where the theory can explain the 
large friction observed in the so called "tacky region" for 
nearly dry contacts, where the capillary bridges pulls the 
rubber into intimate contact with the glass substrate over 
a large fraction of the nominal contact area. 
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